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We investigate the interaction between dark energy and dark matter in the framework of ir- 
reversible thermodynamics of open systems with matter creation/annihilation. We consider dark 
energy and dark matter as an interacting two component (scalar field and "ordinary" dark mat- 
ter) cosmological fluid in a homogeneous spatially flat and isotropic Friedmann-Robertson- Walker 
(FRW) Universe. The thermodynamics of open systems as applied together with the gravitational 
field equations to the two component cosmological fluid leads to a generalisation of the elementary 
dark energy-dark mater interaction theory, in which the decay (creation) pressures are explicitly 
considered as parts of the cosmological fluid stress-energy tensor. Specific models describing coher- 
ently oscillating scalar waves, leading to a high particle production at the beginning of the oscillatory 
period, and models with a constant potential energy scalar field are considered. Furthermore, exact 
and numerical solutions of the gravitational field equations with dark energy-dark matter interaction 
are also obtained. 

PACS numbers: 04.50.Kd, 04.20.Cv, 04.20.Fy 
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I. INTRODUCTION 



The observational data from type la supernovae, ini- 
tially reported in [l|, have generated a large theoreti- 
cal and observational effort for the understanding of the 
observed present accelerated expansion of the Universe. 
Subsequent work on type la supernovae (H, the cosmic 
microwave background and baryon acoustic oscilla- 
tions Q fully support the initial interpretation of the 
observational data that the expansion of the Universe is 
accelerating. The late-time cosmic acceleration is usu- 
ally assumed to be driven by a fluid/field generically 
denoted dark energy Q. Presently very little is known 
about dark energy, namely, its possible composition or its 
structure. Two main scenarios have been proposed to ex- 
plain the nature of the dark energy: a cosmological con- 
stant A 0, or a scalar field, usually called quintessence 
0. The action for gravity and the scalar field 
is S = J[R/16nG-il/2)V"(t)Va(t)-V{^]y^d^x, 
where V{(j)) is the self-interaction potential [3]. 

Another of the central issues in modern astrophysics 
is the dark matter problem (see @ for an extensive re- 
view of the recent results of the search for dark matter) . 
The necessity of considering the existence of dark mat- 
ter at a galactic and extragalactic scale is required by 
two fundamental observational evidences: the behavior 
of the galactic rotation curves, and the mass discrep- 
ancy in clusters of galaxies, respectively. On the galac- 
tic/intergalactic scale the rotation curves of spiral galax- 



ies [SHUj provide compelling evidences pointing towards 
the problems Newtonian gravity and / or standard general 
relativity has to face at these scales. The behavior of the 
galactic rotation curves and of the virial mass of galaxy 
clusters is usually explained by postulating the existence 
of some dark (invisible) matter, distributed in a spherical 
halo around the galaxies. The dark matter is assumed to 
be a cold, pressure-less medium. Many possible candi- 
dates for dark matter have been proposed, the most pop- 
ular ones being the WIMPs (Weakly Interacting Massive 
Particles) (for a review of the particle physics aspects of 
the dark matter see While extremely small, their 

interaction cross sections with normal baryonic matter, 
are expected to be non-zero, so that their direct experi- 
mental detection may be possible. 

In this context, cosmological evolution and dynamics 
are largely dominated by dark energy and dark matter. 
Dark energy has a repulsive effect, driving the Universe 
to accelerate, while dark matter is gravitationally at- 
tractive. In the standard approach to cosmology there 
is no interaction between these two components. Since 
the gravitational effects of the dark energy and of the 
dark matter are opposite (i.e., gravitational repulsion 
versus gravitational attraction) and since dark energy 
is very homogeneously distributed, while dark matter 
clumps around ordinary matter, one expects that any 
dynamic interaction between these two dark components 
of the Universe would be extremely weak, or even neg- 
ligible. However, the possibility of such an interaction 
cannot be excluded a priori and, following some early 
proposals [l^l , presently interacting dark matter-dark en- 
ergy models were extensively investigated in the liter- 
In the standard approach, one may model 
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dark energy as a scalar field with energy density and 
pressure p^. while dark matter is described as a matter 
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fluid with density poM and pressure pdm, satisfying an 
equation of state wum = Pdm/pdm = 0. By assum- 
ing a spatially flat Friedman-Robertson- Walker (FRW) 
background with scale factor, a(f), and allowing for cre- 
ation/annihilation between the dark energy (scalar field) 
and the dark matter fluid at a rate Q, the equations de- 
scribing the variations of the dark energy and dark mat- 
ter densities and pdm can be written as flij 



p^ + iH [l + w^) p^ = 
Pdm + 3H{1 + wdm)pdm = 



(1) 
(2) 



respectively, where H = a/ a is the Hubble function. 
The derivatives with respect to the cosmological time, 
i, will be indicated in the following by an overdot. The 
dark energy and the dark matter create/decay into one 
another via the common creation/annihilation rate ±Q. 
Hence Q describes the interaction between the two dark 
components of the Universe. When Q > 0, dark en- 
ergy is converted to dark matter, while if Q < 0, dark 
matter is converted to dark energy. Since there is no 
fundamental theoretical approach that may specify the 
functional form of the coupling between dark energy 
and dark matter, presently coupling models are neces- 
sarily phenomcnological, although one might view some 
couplings more physical or more natural than others. 
Hence a large number of functional forms for Q have 
been proposed, and investigated in the literature, such 
a,s Q — Pcriti^ + ^)^H{z)Iq{z), where z is the redshift, 
and Iq{z) is an interaction function that depends on the 

redshi ft I l5| . or Q cx pdm4>; and Q cx HpuM, respec- 
tively [Tj"^ 

In writing down Eqs. IJ) and ([2]), one assumes the 
idealised picture of a very quick decay in which the fi- 
nal products reach equilibrium states immediately. How- 
ever, the dark energy/dark matter creation/annihilation 
period may be characterised by complicated nonequilib- 
rium processes, with a highly noncquilibrium distribution 
of the produced particles, subsequently relaxing to an 
equilibrium state. Thermodynamical systems in which 
matter creation occurs fits in the class of open thermo- 
dynamical systems in which the usual adiabatic conserva- 
tion laws are modified thereby including irreversible mat- 
ter creation [l7j . The thermodynamics of open systems 
were first applied to cosmology in [l7j . Explicit inclusion 
of the matter creation in the matter stress-energy tensor 
in the Einstein field equations leads to a three stage cos- 
mological history, starting from an instability of the vac- 
uum. During the first stage the Universe is driven from 
an initial fiuctuation of the vacuum to a de Sitter phase, 
and particle creation occur. The dc Sitter phase docs 
exist during the decay time of its constituents (second 
stage) and ends, after a phase transition into the usual 
FRW Universe. The phenomcnological approach of 
was further discussed and generalised in [1^ through a 
covariant formulation allowing specific entropy variation 
as usually expected for non-equilibrium processes. Cos- 
mological models involving irreversible matter creation 
have been considered in 11911. 



It is the purpose of the present paper to apply the ther- 
modynamics of open systems as developed in [l7j and [isj 
to a cosmological fluid mixture consisting of two compo- 
nents: dark energy, described by a scalar field, and dark 
matter, modeled as an ordinary matter fiuid, in which 
particle decay and production occur. This situation may 
be specific to both early and late stages of cosmologi- 
cal evolution. The thermodynamics of irreversible pro- 
cesses as applied to cosmological models with interacting 
dark energy and dark matter leads to a self-consistent 
description of the dark energy and dark matter particle 
creation/annihilation processes, which in turn determine 
the whole dynamics and future evolution of the Universe. 

The present paper is organised as follows. In Section ITTl 
we present, in some detail, due to its important role, the 
thermodynamical theory of irreversible matter creation 
processes. The theory is applied to a two-component 
cosmological fiuid with interacting dark energy and dark 
matter, and the resulting gravitational field equations are 
written down in Section lHll Particular models, and exact 
and numerical solutions to the field equations are consid- 
ered in Section HVl In Section |V] we discuss and conclude 
our results. Throughout the paper we use a system of 
units so that SttG = c = 1. 



II. THERMODYNAMICS OF IRREVERSIBLE 
COSMOLOGICAL MATTER CREATION 

We consider a cosmological volume element V contain- 
ing particles. For a closed system, N is constant, and 
the corresponding thermodynamic conservation of the in- 
ternal energy E is expressed by the first law of thermo- 
dynamics as [13 



dE^dQ- pdV, 



(3) 



where dQ is the heat received by the system during time 
dt, p is the thermodynamic pressure, and V is any co- 
moving volume. By introducing the energy density p 
defined as p = E /V , the particle number density n given 
by n = N/V , and the heat per unit particle dg, with 
dq ^ dQ/N, Eq. (O becomes 



d { —\ = dq — pd ( — 



(4) 



Equation ^ is also valid for open systems in which N 
is time dependent, N = N{t). Equation (jS]) must be 
modified to explicitly take into account the variation of 
the particle number. 



d (pV) = dQ- pdV + {-]d (nV) 



(5) 



where h ~ p + p is the enthalpy (per unit volume) of 
the fluid. For open systems, adiabatic transformations 
{dQ = 0) are defined by the equation 



d{pV) +pdV - ( ^ ) dinV) = 0. 



(6) 
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Equation (j6]) can be written in a number of equivalent 
forms as 



or, equivalently, 



P = P- P- 



The Einstein field equations 



f 



(7) 



(8) 



(9) 



involve the macroscopic stress-energy tensor T^,y, which, 
in the cosmological case, corresponds to a perfect fluid, 
ft is characterised by a phenomenological energy density 
p and pressure p, and its components are given by 



(10) 



In addition to the Einstein field equations we have the 
Bianchi identities, which lead to Vi/T^ = 0, and to the 
relation 



d{pV) = -pdV. 



(11) 



In the presence of adiabatic irreversible matter cre- 
ation the appropriate analysis must be performed in the 
context of open systems. This involves the inclusion of 
a supplementary creation/annihilation pressure Pc, as we 
may write Eq. ([5]) in a form similar to Eq. (jlip . namely 

m 



dipV)^-ip + Pc)dV, 
where the creation pressure is given by 

d{nV) 



Pc 



dV 



(12) 



(13) 



Creation of matter corresponds to a (negative) supple- 
mentary pressure Pc, which must be considered as part 
of the cosmological pressure p entering into the Einstein 
field equations (decaying of matter leads to a positive 
decay pressure). 



p = p + Pc 



(14) 



The entropy change dS in an open thermodynamic sys- 
tem can be decomposed into an entropy flow d^S^ and 
the entropy creation diS, 



dS = doS + diS, 



(15) 



with diS > 0. To evaluate dS we start from the total 
differential of the entropy, 



Td{sV) = d{pV)+pdV ~ ndinV), 



(16) 



where s = S/V > and pn = h — Ts, p > being the 
chemical potential. In a homogeneous system doS = 0, 
but matter creation contributes to the entropy produc- 
tion. From Eqs. ([15]) and we obtain [l3| 



r^dS 
dt 



T- 



, diS 
~dt 



r^sd{nV)^ 
n dt 



(17) 



To complete the problem we need one more relation be- 
tween the particle number n and V, describing the time 
dynamics of n as a result of matter creation (decay) pro- 
cesses. This relation is usually taken in the form pFzl. [Tsl 



1 d{nV) 
V dt 



(18) 



where ^'(i) is the matter creation (or decay) rate (^'(i) > 
corresponds to particle creation, while 4'(t) < corre- 
sponds to particle decay) . The creation pressure ([T^ de- 
pends on the matter creation (decay) rate, thereby cou- 
pling Eqs. ([T^]) and ([T5)) to each other and, although 
indirectly, both of them with the energy conservation 
law p2p . which is contained in the Einstein field equa- 
tions themselves. The entropy production can also be 
expressed as a function of the matter creation rate as 



S{t) = S{to)ek—'''. 



(19) 



III. COSMOLOGICAL DYNAMICS IN A 
UNIVERSE WITH IRREVERSIBLE DARK 
ENERGY-DARK MATTER INTERACTION 

We shall model the Universe as an open thermodynam- 
ical system, consisting of a two-component (dark energy 
and dark matter) perfect fiuid, with the particle number 
densities denoted by n^, and udm, respectively, cor- 
responds to the "particles" of the scalar field, while tidm 
is the particle number of the dark matter. We denote the 
corresponding energy densities by and pdm, respec- 
tively. The stress-energy tensor of the two-component 
cosmological fluid is given by 

T; = T^.^)'' + Tl'^^'^ = pu^u^ - p 5';, (20) 

where m'' = dx'^/ds is the four- velocity, and 



P = P4> + Pdm, 



P 



Pdm- 



(21) 



The energy density and pressure of the dark energy 
are given by p^ = (/>^/2 -I- U{4>) and p^ — (j)'^/2 ~ U{(j)), 
respectively, where U{(j)) is the self-interaction potential. 
We suppose that neither the particle numbers nor the 
stress-energy of the components are separately conserved, 
that is, particle inter-conversion and exchange of energy 
and momentum between the two components are admit- 
ted. The cosmological fluid mixture is characterised by 
a total energy density p — p^ + Pdm, total thermody- 
namic pressure p = Pc/, + Pdm and a total particle num- 
ber n = + riDM- We consider that the geometry of 
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the spacetime is described by the flat FRW line element, 
given by 



ds^ = dt^ - a^{t) {dx^ + dy^ + dz^) 



(22) 



We shall assume that the particle number densities 
and nuM of each component of the fluid obey the follow- 
ing balance laws, 



(23) 
(24) 



respectively, where Fi 7^ and V2 ^ ^ are arbitrary 
functions. Equations (|23p and p4)) describe the decay of 
the dark energy 0-particles, and the creation of the dark 
matter particles, with a scalar field decay rate and a dark 
matter creation rate ^{t) (x p^. Thus, the dynamics 
of 0-particle decay and the creation of dark matter is 
governed in the present model by the scalar field via its 
energy density. From Eqs. (|23|) and ()24p it follows that 
the total particle number n obeys the balance equation 



3Hn,= (Fa-Fi), 



(25) 



Hence in the case of an interacting dark matter and dark 
energy the total particle number conservation occurs only 
in very special cases, and therefore we shall suppose that 
generally Fi ^ F2. 

In the framework of the thermodynamics of irreversible 
processes particle creation and decay gives rise to a decay 
and a creation thermodynamic pressure, given by 



(0) 



Ti jPci, + P4>) P<i, 



(26) 



and 



P. 



(DM) _ r2 {pDM +Pdm) P4, 



iHuDM 

respectively, while the total creation pressure becomes 



(27) 



P, 



(total) 



(0) 



P: 



3H 



P. 



(DM) 



Ti {p^ + P4,) T2 {pDM + Pdm) 



(28) 



Using the results obtained above the complete Einstein 
gravitational field equations describing the dynamics of 
a flat FRW spacetime filled with a mixture of interact- 
ing dark matter (scalar field) and dark matter can be 
expressed in the form 



2H 



m — P4, + Pdm, 

= - Pdm 



3H 



^2 {PDM +PDm) 



ricf, 



and 



fi^ + 3Hn^ 

flDM + iHuDM 



(29) 



(30) 

(31) 
(32) 



where poM = Pdm {noM) and pdm = Pdm {jidm), and 
recalling that the energy density and pressure of dark 
energy is given by 



P0 = y + C/(0), P4> = \-U{<j>). (33) 

As applied to each component of the cosmological fiuid, 
Eq. (O , the second law of thermodynamics for open sys- 
tems provides the following relationships 



+ m {p^ +p^) + -Ll 



+ 1 



0, 



(34) 



: QUI- I A ^2{PDM +PDm) P<p 

PDM + aH [pdm + Pdm) = , (35) 



udm 



respectively. 

Equation ([M)). which describes the dynamics of the 
dark energy during its interaction with dark matter, can 
be written in an equivalent form as 



'(j) + 3i70 + F (0, (/),[/)</) + U'{(j)) = 



(36) 



where we have denoted F y(/),(f>,Uj ~ Vxp^jn^. There- 
fore in the framework of the thermodynamics of ir- 
reversible processes a friction term in the scalar field 
Eq. p6| arises naturally, and in a general form, as a 
direct consequence of the second law of thermodynamics 
as applied to an open system. 

Adding Eqs. and the evolution of the total 
energy density p = p^ + Pdm of the cosmological fluid is 
governed by the equation 



P + 3H (p + p^+Pdm) 
^2 (pdai +Pdm) 



El 



(37) 



noM 

For the entropy of the newly created matter we obtain 



SDMit) = Sdm (io)exp I / F2 

/to noAi 



dt 



(38) 



Irreversible matter particle creation is an adiabatic 
process, the produced entropy being entirely due to the 
increase in the number of fluid particles, there being no 
increase in the entropy per particle due to dissipative 
processes. 



IV. IRREVERSIBLE DARK ENERGY-DARK 
MATTER INTERACTION MODELS 

In the present Section we consider, within the frame- 
work of irreversible thermodynamics with matter cre- 
ation/annihilation, a number of speciflc cosmological 
models with dark energy-dark matter interaction. As a 
first case we consider the situation in which the density of 
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the dark matter is much smaller than the energy density 
of the scalar field. This case corresponds to an Universe 
dominated by the dark energy component, assumed to 
be represented by a coherent wave of 0-particles. In this 
case the kinetic term dominates in the total energy of 
the scalar field. As a second case we consider a poten- 
tial energy dominated scalar field. The general dynamics 
of the cosmological model with interacting dark energy 
and dark matter is also considered, and the cosmological 
evolution equations are studied numerically. 



A. Coherent Scalar Waves-Dark matter interaction 

We shall consider in the following that the energy den- 
sity and particle number of the newly created dark mat- 
ter is much smaller than the energy density and parti- 
cle number of the corresponding scalar field fluid compo- 
nent, that is the relations poM ^ P<t>, noM ^ n^j,, and 
Pdm Pct> hold. In this case the Universe is dominated 
by the scalar field energy density, and its evolution is not 
influenced by the matter content. We shall work through- 
out with finite values of the fluid quantities at < — to- The 
coupling between scalar field and dark matter is realised 
only by means of the balance equation of ordinary matter 
via the scalar field energy density, and the basic equations 
describing the dynamics of a flat FRW scalar fleld filled 
space-time interacting with a dark matter component are 
given by 



3Hn^ 



and 



P4> 



(39) 
(40) 

(41) 



(42) 
(43) 



respectively. 

A homogeneous scalar field oscillating with frequency 
can be considered as a coherent wave of "particles" 
with zero momenta, and with a particle number density 
given by [2^ 



— constant. 



(44) 



In other words, oscillators of the same frequency 
oscillating coherently with the same phase can be de- 
scribed as a single homogeneous wave 4>{t). Insertion of 
the energy density of the scalar field given by Eq. p4|) in 
Eq. (pij) leads to the condition 



0, 



(45) 



or, equivalently. 



U 



(*)4- 



(46) 



Therefore, a homogeneous oscillating scalar field is de- 
scribed in the present model by a Barrow-Saich type po- 
tential, with potential energy of the scalar field propor- 
tional to the kinetic one [2l| . The energy density of the 
scalar field becomes p^ = (p^, and this relation, obtained 
naturally in the framework of the present formalism is 

very similar to the equation p^ = {'^^) obtained by re- 
placing by its average value per cycle [2^ . 

In this case the equations describing the dynamics of 
the FRW type spacetime filled by the decaying oscillating 
homogeneous scalar field in the presence of dark matter 
creation become 



2H + 3H'^ = -Tim^H, 



nuM + 3Hn 



DM 



(47) 



(48) 



(49) 



respectively. In the following we will assume, for simplic- 
ity, that Fi and F2 arc constants. By introducing a set 
of dimensionless variables r. h and 6om by means of the 
transformations 



t ■ 



H 



Fim^ 



TT-DM 



2FiF2TOrf 



Fim^ 3 
Equations (gS]) and dH]) take the form 

-h{h + l), 



dh 
d7 



d0 



DM 



dr 



and yield the following general solutions 

1 



h{T) 



and 



TDM 



(e^° - 1)^ 



^DMQ 



+ 1/2 e2(^-^«) - 1/2 



'DM, 

(50) 
(51) 
(52) 

(53) 
(54) 



respectively, where we have denoted 9dmo = 0dm {tq)- 
The evolution of the scale factor is given by 



a(T) = ao ( 1 



(55) 



The deceleration parameter q = d{l / H) / dt — 1 is given 
by 



(56) 



6 



For {t - to 



i.e., at the start of the osciUatory 



-TUl yi^ — LQ) ^ X , i.e., clt tile SLUit Ui tlie U&CiiicltUiy 

period corresponding to the dark matter production, the 
approximate solution of the field equations is given by 



aoT 



2/3 



To 



^dmqTq 



(57) 



(58) 



This phase corresponds to an Einstein-de Sitter expan- 
sion, with decaying dark energy, and dark matter cre- 
ation. 

During the initial oscillating period of the scalar field 
dominated FRW Universe there is a rapid increase of 
its dark matter content. The particle number density 
increases during a time interval 



At = tr 



in 



1 



2r, 



to, 



and reaches a maximum value given by 



(max) 



1 - 361 



3At ^o(l-3^^I5Mo^o/4^2)■ 



(59) 



(60) 



Assuming that the scalar field oscillations decay into 
rclativistic dark matter, the energy density and temper- 
ature of the dark matter component of the cosniological 



fiuid is given by poM 



DM 



'i' ^, where 7 = 2 



corresponds to a stiff (Zcldovich) fiuid obeying an equa- 
tion of state of the form pdm ~ Pdm, and 7 = 4/3 
corresponds to a rclativistic, radiation like fiuid. 

The entropy produced during dark matter particle cre- 
ation can be easily obtained from Eq. (|38p and is given, 
in first approximation, by 



T -To 



Sdmo OdmoTq 



1. 



(61) 



Therefore, for small times, there is a linear increase of 
the entropy of the dark energy dominated fiat spacetime. 



Constant scalar field potential dark energy and 
darli matter interaction 



As a second example in the study of the irreversible 
interaction between dark energy (a scalar field) and dark 
matter we consider the case in which the scalar field po- 
tential may be approximated, at least for a certain time 
interval, as a positive constant, U (0) = A = constant > 
0. Therefore, the energy density and the pressure of the 
scalar field can be written as 



2^ 



■A, 



(62) 



Moreover, we assume that the change in the scalar field 
energy density due to the cosmological expansion can be 



neglected in the scalar field evolution equations, that is, 
in the energy density and particle number equations of 
the scalar field, given by Eqs. (|3T|) and (|34|) . the term con- 
taining 3i? can be neglected. In this approximation the 
main contribution to the temporal dynamics of the scalar 
field is due to its decay into dark matter particles, and not 
to the cosmological expansion. On the other hand, the 
feedback of the newly created dark matter particles on 
the cosmological dynamics cannot be neglected. Within 
this approximation from Eq. pip , describing the scalar 
field particles decay, we obtain first 



i 1 



(63) 



which gives the scalar field particle number density as a 
function of the energy of the scalar field. By substituting 
this expression of into Eq. gives 



ns + El — fiA, = 0. 



Taking into account that 



= 2 (p^ - A) 



T7 



A 



(64) 



(65) 



we obtain the equation describing the dynamics of the 
scalar field particles as 

n^n^ - 2ri^ - ^TiKh^ = 0. (66) 

By denoting = m, fi^ = udu/dn^, Eq. (|66| becomes 

du 



providing 



n^— =2(u + riA). 
dn^ 



n0 = Nin^ - EiA, 



(67) 



(68) 



where A^i is an arbitrary integration constant, which 
can be determined from the initial condition (to) = 
-Fip^ (to) = -Eip^o and {to) = n^o as 



Fi (A - p^o) 



(69) 



Therefore the general solution of Eq. (j66p is given by 
/r\A 1 + iV2e2"* 



nff,{t) 



Ni 1 - iVse 



2at ' 



where a = vTiAJVi , and 



n^o + x/FiA/A^i 



(70) 



(71) 



By assuming that the newly created dark matter parti- 
cles are pressureless, pdm = 0, and the energy density of 
the dark matter component is Pdm = 'mDM'n-DM, where 



iriDM is the mass of the dark matter particle, then the 
evolution equation for the dark matter density can be 
written as 

^ 1 



-2 "0- 



(72) 



By neglecting the effects of the expansion of the Universe 
the dark matter energy increases as 



0.052 - 



PDAlit) 



To A 



r2(A- 



riA A - 
Ni n2 



riA 



(73) 



In the general case, in which the expansion of the Uni- 
verse is also taken into account, the system of equations 
describing the irreversible dark energy-dark matter inter- 
action are given by 



(74) 



FIG. 1: Time variation of the scale factor of the interacting 
dark energy-dark matter filled Universe, for different values 
of the parameter A: A = 0.001 (solid curve), AO. 005 (dotted 
curve), A = 0.009 (dashed curve), and A = 0.012 (long dashed 
curve). The initial conditions used for the numerical integra- 
tion of the cosmological evolution equations are a(0) = 0.05, 
N^{Q) = 10, r0(O) = 1.5, and rDM{Q) = 0.01, respectively. 
For the mass moM of the dark matter particle we have as- 
sumed the value m = 1 MeV. 



(75) 



p, + 6H ip, - A) + ^^^(P^^-^^P^ = 0, (76) 



Pdm + SHpDM — TimDMP<t>-, 



(77) 



where, for simplicity, we have assumed Fi = r2 = 
constant. By introducing a set of dimensionless variables 
(r^, A^0,roA/,T), defined as 



p^ ^ Ar^, ^ AiV^, 
Pdm = AmDAirDM, t = T/Ti, 



(78) 
(79) 



the field equations Eqs. (|74p - ([77)) can be written in a 
dimensionless form as 



1 da . 

-3- = A^r0 + mDAirDM, 
adr ^ 



(80) 



dr 



+ 3X\/r^ + mDAirDM = -r^, (81) 



'^''^ , c\ / \ f T\ , 2r0(r0-l) 

+ bX^Tft, + niDMrDM [Tcf, - 1) H — = 0, 



dr 



(82) 



dr 



DM 



dr 



+ 3A^r0 + niDAirDAirDAi = r^, (83) 



respectively, where 



(84) 




FIG. 2: Time variation of the scalar field particle number A^,^ 
in the interacting dark energy-dark matter filled Universe, for 
different values of the parameter A: A = 0.001 (solid curve), 
AO. 005 (dotted curve), A — 0.009 (dashed curve), and A = 
0.012 (long dashed curve). The initial conditions used for the 
numerical integration of the cosmological evolution equations 
are a(0) = 0.05, iV0(O) = 10, r^(0) = 1.5, and roAiiO) = 0.01, 
respectively. For the mass nioM of the dark matter particle 
we have assumed the value m = 1 MeV. 



The dynamics of the interacting dark energy-dark matter 
system is determined by two control parameters. A, and 
the mass mDM of the dark matter particle. The time 
variations of the scale factor a, of the scalar field particle 
number N^, of the scalar field energy r^, as well as the 
dark matter energy density tdm are represented, for dif- 
ferent values of the parameter A, and for ttldai = 1, in 
Figs.IIEl 
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FIG. 3; Time variation of the dimensionless scalar field en- 
ergy in the interacting dark energy-dark matter filled Uni- 
verse, for different values of the parameter A: A = 0.001 (solid 
curve), AO. 005 (dotted curve), A = 0.009 (dashed curve), and 
A = 0.012 (long dashed curve). The initial conditions used 
for the numerical integration of the cosmological evolution 
equations are a(0) — 0.05, N^{0) — 10, '■^(O) — 1.5, and 
rDAiiO) = 0.01, respectively. For the mass vtldm of the dark 
matter particle we have assumed the value m = 1 MeV. 




1 2 3 4 5 6 7 



FIG. 4: Time variation of the dimensionless dark matter en- 
ergy TDM in the interacting dark energy-dark matter filled 
Universe, for different values of the parameter A: A = 0.001 
(solid curve), A0.005 (dotted curve), A = 0.009 (dashed 
curve), and A = 0.012 (long dashed curve). The initial con- 
ditions used for the numerical integration of the cosmological 
evolution equations are a(0) = 0.05, N^{0) = 10, r^{0) = 1.5, 
and r£)A/(0) = 0.01, respectively. For the mass moM of the 
dark matter particle we have assumed the value m — 1 MeV. 

As one can sec from Fig.[Tl the interacting dark energy- 
dark matter filled Universe is an expansionary state, with 
the rate of the expansion, and the scale factor evolution, 
strongly dependent of the dimensionless parameter A, 
which depends on the ratio of the (constant) scalar field 
potential A, and of the scalar field decay rate Fi. Accel- 
erated expansion can also be obtained in the framework 
of the present model, the creation pressure correspond- 



ing to the irreversible decay of the scalar field, and the 
matter creation, can drive the Universe into a de Sitter 
type phase. The scalar field particle number decays dur- 
ing the cosmological evolution, as shown in Fig. [21 The 
decay rate strongly depends on the numerical value of A. 
The dimensionless energy of the scalar field r^, shown 
in Fig. [21 tends in the large time limit to the value 1, 
corresponding to — A, and to a de Sitter expansion. 
This shows that the decay of the scalar field is deter- 
mined and controlled by the kinetic energy term of the 
field <p^/2, which is the source of the dark matter parti- 
cles creation. When the energy and the pressure of the 
scalar field are dominated by the scalar field potential A, 
Pcf, = —pcji = A, then + P4, ~ 0, and from Eq. ([M)) it 
follows that p0 — constant, and the scalar field energy 
cannot be converted any more into other type of parti- 
cles. The energy density of the dark matter particles, 
presented in Fig. [31 increases in time due to the decay of 
the (/)-particles. 



V. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have shown that the thermo- 
dynamics of the open systems is a valuable tool for de- 
scribing the interacting dark energy - dark matter phases 
of the general relativistic cosmological models. As ap- 
plied to a two component (scalar field and dark matter) 
cosmological model, the thermodynamics of irreversible 
processes provides a generalization of the elementary the- 
ory of dark energy-dark matter interaction, which envis- 
ages that during the dark energy dominated phase of 
the expansion of the Universe, when the expansion slows 
down, the energy stored in the zero mode oscillations of 
the scalar field transforms into particles via single parti- 
cle decay. Thus the model presented in this paper gives 
a rigorous thermodynamical foundation, and a natural 
generalisation, to the theory of dark energy-dark mat- 
ter interaction. Particle decay (creation) gives rise to a 
supplementary decay (creation) pressure which has to be 
included as a distinct part in the stress-energy tensor of 
the cosmological mixture. We have considered only par- 
ticle creation in a scalar field (dark energy) dominated 
Universes, leading to a model in which particle produc- 
tion rate is extremely high at the beginning of the oscil- 
latory period, and afterwards it tends to zero, when the 
scalar field energy and pressure become dominated by 
the scalar field potential, assumed to be a constant, and 
when the kinetic energy of the scalar field becomes negli- 
gibly small. Such a dark matter particle production has 
very important implications on the dynamics and evolu- 
tion of the Universe. The details of the dark energy-dark 
matter interaction mechanism depend on the parameters 
of the particle physics models involved to describe the 
newly created particles. One of the parameters is the 
dark matter particle mass, which is a key parameter in 
the description of the scalar field-dark matter process. 
Unfortunately presently there is no definite answer giv- 
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ing the value of the mass of the dark matter particle. 

Particle creation can be related to what is called the 
arrow of time: something that provides a direction to 
time, and distinguishes the future from the past. There 
are two different arrows of time: the thermodynamical 
arrow of time, the direction in which entropy increases, 
and the cosmological arrow of time, the direction in which 
the Universe is expanding. Particle creation introduces 
asymmetry in the evolution of the Universe, and enables 
us to assign a thermodynamical arrow of time, which 
agrees, in our model, with the cosmological one. This 
coincidence is determined in a natural way by the decay 
of the scalar field, due to the presence of a friction force 
between field and matter. 

In our present approach we have neglected the back- 



reaction of the newly created particles on the dynamics 
of the Universe, and the effect of the form of the scalar 
field potential U {(/)) have not been completely envisaged. 
These aspects of the thermodynamic theory of the dark 
energy dark matter interaction will be the subject of a 
future work. 
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